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ABSTRACT 

More than two dozen short-period Jupiter-mass gas giant planets have been 
discovered around nearby solar-type stars in recent years, several of which un- 
dergo transits, making them ideal for the detection and characterization of their 
atmospheres. Here we adopt a three-dimensional radiative hydrodynamical nu- 
merical scheme to simulate atmospheric circulation on close-in gas giant planets. 
In contrast to the conventional GCM and shallow water algorithms, this method 
does not assume quasi hydrostatic equilibrium and it approximates radiation 
transfer from optically thin to thick regions with flux-limited diffusion. In the 
first paper of this series, we consider synchronously-spinning gas giants. We show 
that a full three-dimensional treatment, coupled with rotationally modified flows 
and an accurate treatment of radiation, yields a clear temperature transition at 
the terminator. Based on a series of numerical simulations with varying opacities, 
we show that the night-side temperature is a strong indicator of the opacity of the 
planetary atmosphere. Planetary atmospheres that maintain large, interstellar 
opacities will exhibit large day-night temperature differences, while planets with 
reduced atmospheric opacities due to extensive grain growth and sedimentation 
will exhibit much more uniform temperatures throughout their photosphere's. In 
addition to numerical results, we present a four-zone analytic approximation to 
explain this dependence. 



1. Introduction 

Among the over 200 extra-solar planets discovered through a variety of observational 
techniques including Doppler shifted spectral lines, transit searches, and gravitational lens- 
ing, approximately 9% have orbital periods less then 4 days. One of the more puzzling and 



unexpected results to come from planet detections, these planets, known as hot- Jupiter's 
or Pegasus planets, receive approximately 10 4 times more energy from stellar irradiation 
than from internal heating. This energy input greatly modifies the thermal structure and 
dynamics of the atmosphere. Given the short orbital period of these planets, tidal forces also 
play a stron g role, quickly circularizing their orbits and synchro nizing; their spin and orbital 
frequencies (IGoldreich Sz Soterl Il966l ; iDobbs-Dixon et al.l 120041 ) . Consequently, one side of 
the planet perpetually faces the host star, reaching equilibrium temperatures of over 1000K, 
while night-side temperatures are determined by the ability of the planet to redistribute 
energy. In the absence of heat redistribution, the night-side of the planet would remain 
at « 100K after a few Gyr of cooling and Kelvin-Helmholtz contraction. For comparison, 
Jupiter has a rotational period of 11 hours and its internal luminosity is approximately equal 
to the energy it receives from the Sun. Its photospheric temperature is quite uniform in both 
longitude and latitude. 

The small semi-major axis of hot- Jupiter's provide several new m ethods for exploring 
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1. Detection of a primary eclipse, as the planet 



passes between the star and the Earth, yields measurements of planetary radii and masses 
that can be used to infer structural parameters of the planets. The secondary eclipse, caused 
by the passage of the planet behind its host star, provides a direct measurement of the 
day-side temperature of the objects by measuring the decrement in infrared flux. As of the 
writing of this paper, 14 planets have been detected via the primary transit method, while 
3 of those have also been detected via their secondary eclipse. The day-side temperatures of 
the objects f ound from secondary eclipse measurements are found to be T = 1130 ± 150K fo r 
HD209458b foeming et alil2005bh . T = 1060 ± 50K for Tr ES-1 Jcharbonneau et aljlioosl l. 
and T = 1117 ± 42K for HD 189733b foeming et alil2006h . Measurements of HD209458b 



and HD189733b were done using the Spitzer Space Telescopes 24/wn and 16/xm band-passes 
respectively. The estimate of temperature is then sensitive to both the ratio of stellar and 
planetary radii and stellar temperature. TrES-1 was observed in both the 4.5/xm and 8/im 
bands, and temperatures were derived by assuming the planet emits as a blackbody. In 
addition, recent non-transiting observations of v— Andromeda b with the S pitzer Space Tele- 



scope , indicate a substantial orbital phase dependence in the flux at 24/im ( [Harrington et al. 
2006h . 



A number of groups have used one-dimensional, pla ne-parallel radiative t ransfer models 



to ca l culate the emergen t spectra of these hot - Jupiter's (IFortnev et a. 
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transfer, radiative equilibrium, and hydrostatic balance to determine the pressure, temper- 
ature, and spectral energy distribution as a function of radial position in the atmosphere. 
The input parameters include the temperatures at the bottom of the atmosphere and the 
incident stellar energy spectrum at the top. The temperature at the bottom is determined 
primarily by the internal luminosity, and ideally it should be taken from one-dimensional 
evolutionary and structural models. However, given the dominate influence of the incident 
energy near the surface, the choice of the bottom does not noticeably influence the emergent 
spectra. 

In hot- Jupiter's, the intense irradiation of the day-side drives strong thermal winds to- 
ward the night-side. The resulting temperature in the upper atmosphere depends on the 
ability of these winds to redistribute the stellar irradiation, and in general should be a func- 
tion of both longitude and latitude. In the absence of dynamical models, this redistribution 
of incident energy at the upper layers is usually set to be some fraction of the incident irra- 
diation, to represent the degree of energy re-distribution. The value of this parameter is a 
major uncertainty in these radiative models, and authors have computed a number of cases 
with varying degrees of re-distribution. In addition, such parameterization neglects the role 
of advection and radiative transfer within lower levels of the atmosphere, which are also 
important in determining the final pressure-temperature profiles. 

In an attempt to address the dynamical redistribution of incident stellar energy within 
the atmosphere, there have been a number of dynamical models of the atmospheres of hot- 
Jupiter's. These models utilize a variety of methods employ i ng various simplifications, 
inclu ding; solving the primitive eq uations (IShowman fc Guillotl 2002 ; ICooper fc Showman 



2005 ). the shallow water equations fjCho et al 



2003 



Menou et al. 



2003 



Langton fc Laughlin 



20071 ). the equivale nt-barotropic equatio ns flCho et al.ll2006l ). and two-dimensional hydrody- 
namical equations (IBurkert et al.l 120051 ). These models predict a wide range of behavior. 



More details on these mod els are presented in jjH in comparison with the model presented 
in this paper. One group (IFortney et al.ll2006a ) has attempted to couple the spectral and 
dynamical mo dels by utilizing pres s ure a nd temperature profiles derived directly from the 
simulations of ICooper fc Showman! (120051 ). Unlike previous models, these atmospheres are 
not iterated to achieve radiative-convective equilibrium. 

In this paper, we present the results obtained from a three- dimensional, hydrodynamic 
simulation based on flux-limited radiative transfer models of hot- Jupiter's for a variety of 
rotation rates and opacities. In §0 we present the basic equations, numerical methods, and 
initial conditions. In £j3]we present our results and analysis for both rotating and non-rotating 
flows. We also study the effects of changing opacity on the dynamics and heat distribution. 
In §4] we include a detailed comparison with previous dynamical models in an attempt to 
highlight the consequences of making certain simplifying assumptions. We conclude in §5] 
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with a discussion. 



2. Numerical Method 

2.1. Flux-Limited Radiative Hydrodynamical Model 

We cast our numerical model in spherical coordinates (r, <p, 9), where <fi is the azimuthal 
angle (or longitude) and 9 is the meridional angle (or latitude) measured from the equator. 
Including both the Coriolis (2f2 x u) and centrifugal forces (Ox (ft x r)), the equations of 
motion for the fluid can be written as 

^ + V-(pu) = (1) 



+ (u ■ V) u = — VP + -V$ - 2Q x u - Ox (n x r) (2) 

ot p p 

The rotation frequency is given by f2, and the gravitational potential, $ = — , varies 

only in the radial direction. We neglect explicit viscosity, but some degree of numerical 
viscosity is inevitable. We also include the curvature terms in (u • V) u. We solve equations 
([1]) and (|2J) on a stagged grid, where scalars are defined in the center of cells and vectors on 
cell boundaries. This method yields second-order spatial accuracy. Given the decreasing grid 
size near the pole, our computational domain is limited to 9 = ±70°. Although excluding 
this region neglects an avenue for energy re-distribution, even a modest amount of rotation 
(approximately 2 to 4days for the hot- Jupiter's) will cause the dominate flow patterns to 
be concentrated near the equator. We simulate the entire azimuth of the planet, instituting 
periodic boundary conditions at <p = and 2tt. The radial extent of the domain extends 
from 7.95 x 10 9 cm to 8.65 x 10 9 cm, corresponding to 1.06 to \.2Rj. The pressure scale height 
on the day-side is approximately 340 km. 

Our numerical radiative transfer scheme is capable of following the temperature and 
radiation energy independently, linking them with a given heating/cooling function. How- 
ever, this degree of sophistication is not necessary for calculations in which the gas is close 
to thermal equilibrium. Instead, we use a one-fluid approximation where radiation energy 
density is a simple function of temperature, E = aT A . For the temperatures considered here, 
the radiative energy density is much lower than the thermal energy density. The internal 
and radiation energy equations reduce to 



<9e / ,_x 
- + (u.V)e 



-PVu-V-F (3) 
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where e = c v pT is the internal energy density, T is the temperature, c v is the specific heat, 
and F is the radiative flux. 

To calculate the radiat i ve flux we use the flux limited radiation transfer approximation 
of iLevermore fc Pomraningi (1198 IT ). 



F = -X—VE, 



(4) 



where A is a non-constant flux limiter which prescribes the relationship bet ween the radiative 

flux an d the radiative energy gradient. We use the flux limiter developed by ILevermore fc Pomraning 
<jl98lh . given by 

A=^4±^. (5) 
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where 



R 



3R + R 2 ' 
VE\ 



1 



(6) 



pn E 

compares the scale height of the radiation energy to its mean free path. This overall pro- 
cedure allows for an accurate treatment of both optically thick and thin regions; in the 
optically thick diffusion limit R — > 0, A — > 1/3, and F — > — c\7E/3pn, while in the opti- 
cally thin streaming limit R — > oo, A — > 1/R, and F — > cE. While the flux limiter A is an 
approximation, it is quite accurate in both the optically thick and thin limits. 

In order to follow the evolution on a long time-scales, the radiative portion of the energy 
equation is solved implicitly. For quasi-static radiative conditions such as those considered 
here, the equations can be advanced much more rapidly than if they were restricted by a 
radiative time-step. We use the successive over-relaxation method (SOR) to solve the V • F 
portion of equation Q, alternately updating even and odd grid cells. The limiting factor 
for the numerical time-step then becomes the Courant condition. 



The advection scheme, described i nlHawlev et al.l (119841 ) and lKley fc Henslerl (119871 ). is 
an extension of the first-order van Leer (Ivan Leerlll9771 ) method known as the 'mono-scheme'. 
It employs operator-splitting, where the finite difference equations are split into parts, which 
are then solved separately always using the latest values of the variables. The scheme yields 
semi-second order temporal accuracy and allows for the accurate resolution of discontinuities 
in the fluid flow with limited diffusion. 

Finally, we use an ideal gas equation of state for the pressure, p = (RcpT)/ p, with 
specific heat c v = Rc/(/i(7 — 1)) where is the gas constant, 7 = §, and the mean molec- 
ular weight is fixed at p — 2.3. Although the temperatures in the hydrodynamic models do 
become hot enough to dissociate hydrogen molecules in some locations within the planet, 
the region which we are most concerned with is well described by a constant molecular 
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weight. Radiative opaci ties are found using the t ables of iPollack et al.l (119851 ) for lower tem- 
peratures coupled with lAlexander fc Ferguson! (119941 ) for higher temperatures. These are 
Rosseland mean opacities and include the effects of atomic, molecular, and solid particulate 
absorbers and scatters. The opacity is one of the largest uncertainties. The effect of compo- 
sition, clouds, settling, wavelength dependence, grains' condensation, sublimation, collisional 
growth, and sedimentation are but a few of the parameters that alter the magnitude of the 
opacity. To address this uncertainty, we explore the effect of varying opacity in §3.31 

Our initial density and temperature profiles are taken from a one-dimensional pl a netar y 
evolutionary code. For details on the model, see model Bl of iBodenheimer et al.l (120011 ). 
Use of this model implies that the simulations are initially spherically symmetric. Although 
unrealistic, we allow simulations to run for sufficient time to relax into equilibrium config- 
urations; the details of the initial state are lost. Figure ([1]) shows the initial temperature 
and density profiles from the one-dimensional model near the top of the planet. The model 
follows a 0.63M/ planet for 4.5Gyr. The upper boundary was held at a fixed te mperature 



of 1200K to simulate the effects of irradiation from the central star. As noted in Cho et al. 



(120061 ). a quiescent initial start neglects the effects of pre-established small-scale structures 
such as eddies and jets. Given that we lack the resources to simultaneously study these small 
scale effects and full three dimensional effects, we assume that the dynamics will be domi- 
nated by the large scale anisotropic heating imposed by the stellar irradiation. In contrast to 
the shallow-water approximation, specific vorticity can be generated in our full 3D radiative 
hydrodynamic simulations. Therefore, baroclinic instabilities can be excited spontaneously 
and can lead to the generation of structure down to the resolution length scales. 



To represent the impinging radiation, we set the temperature at the upper boundary to 



T 



max 



1200K (cos (0) cos (fl)) 1/4 , 100K 



approximating the stellar radiation field with a 
maximum at the sub-stellar point, and a night-side held at 100K. The night-side tempera- 
ture is chosen to be consistent with the photospheric temperature derived in non-irradiated 
planetary evolution models. This imposed temperature distribution will increase the scale 
height of the atmosphere on the day-side, while cooling and contracting it on the night- 
side. An interesting result is a planet that is no longer spherical; the scale-height on the 
day-side is somewhat larger than on the night-side. Given the computational difficulties of 
modeling the extremely low density regions in the upper atmosphere, we impose a movable 
upper boundary at the location where p < 10 _9 g/cm 3 . At the bottom boundary we specify 
the temperature flux taken from our initial one-dimensional models. As mentioned above, 
the energy input from the star overwhelms the planets intrinsic luminosity in the upper 
atmosphere, so the choice of this boundary condition is not critical. We have run several 
simulations with fixed temperature to verify the validity of this approximation, and found 
little difference in the resulting energy distribution or dynamics. 
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Fig. 1. — Spherically s ymmetric temperatu r e and density profiles from the one- dimensional 
evolutionary models of iBodenheimer et al.l (120011 ) used as an initial condition for the simu- 
lations presented here. The models were run for 4.5Gyr, during which time the temperature 
at the upper boundary was held at 1200K to simulate the irradiation from the central star. 
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3. Results of Numerical Simulations 

3.1. The Non-Rotating Model 

As a reference, we first present the results from an idealized and artificial non-rotating 
model in which the gas giant is subject to one-sided stellar irradiation. The azimuthal 
pressure gradient associated with the imposed temperature contrast drives strong winds 
toward the night-side of the planet. The flow pattern is symmetric with respect to the planet- 
star line, and fluid flows to the night-side around both sides of the planet. As the material 
moves to the night-side, it achieves maximum velocities of ~ 3km/s near the terminators 
(</> = f and (j) = ^f)- This velocity corresponds to Mach numbers of up to 2.7. A substantial 
portion of energy is released in these shocked regions. By the time the two symmetric flows 
converge on the night-side they have both undergone substantial cooling and sink radially 
inward, initiating a return flow at depth. In Figure ([2]) we show the temperature distribution 
in the equatorial plane. Despite the high winds, a clear day/night difference persists. Also 
evident is the cool region at depth, which is due to the confluence of the negative temperature 
gradient from the interior, the positive temperature gradient near the surface, and the cool 
flow returning to the day-side, completing the two, approximately symmetric, azimuthal 
convection cells. 

Figure ([3]) shows the temperature and velocity at the photosphere (r = |) of the planet. 
Given the changing scale height near the surface of the planet, these plots do not represent 
a constant radial surface. Moving from the day to night-side, a clear increase in the velocity 
near the terminator is evident, followed by a drop at the convergent point. It should also 
be noted that the convergence point is quite dynamical, oscillating in both longitude and 
latitude. However, despite these oscillations, two distinct convective cells remain. Night-side 
temperatures range from 300 — 550K, with an average night-side temperature (| < < ^) 
of 414K. Also evident in Figures ([2]) and ([3]) is an area of increased temperature near 180°. 
As the two converging flows meet on the night-side, compressional heating drives up the 
temperatures, and thus the local scale height. The resulting increase is small, but can be 
seen both at depth and near the photosphere. 

3.2. A Planetary Model with 3-Day Rotation 

It is widely believed that tidal forces within the atmospheres and the envelope of hot- 
Jupiter's drive them to tidally locked spin configurations on time-scales much shorter than 
the main sequence life span of their host stars. In this synchronous state, the planets' spin 
frequency equals their orbital frequency. With this assumption, hot- Jupiter's spin with pe- 
riods on the order of 3— days. In comparison to our own giant planets, this is a relatively 
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-1.0 -0.3 0.0 0.S 1.0 



Fig. 2. — The temperature distribution at the equator of a non-rotating hot- Jupiter. The 
inner boundary is assumed to be spherically symmetric, with an outward energy flux fixed 
from the initial one-dimensional model. The outer black area represents regions with p < 
10~ 9 g/cm 3 , outside our movable boundary. 
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6.3e+02 5.6a+04 1.le+05 1.7e+05 2.2e+05 2.8e+05 3.3e+05 



Fig. 3. — The temperature (upper panel) and velocity (lower panel) at the photosphere of a 
non- rotating planet. A strong temperature gradient between the day-side (0 = 0°) and the 
night-side (0 = 180°) is evident. The regions of largest VT correspond to fluid motions of 
> 3km/ s. 
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slow spin rate. Nevertheless, the associated Coriolis force significantly alters the resulting 
flow dynamics and may have implications concerning the ability of the planet to fully syn- 
chronize its spin. Another crucial consideration that came to light during this study was the 
effects associated with the initialization of the rotation. The one-dimensional initial models 
described in Section (12.11) were non-rotating models, and the rate that we chose to turn on 
the rotation had observable effects. We will explore this effect in a subsequent paper. 



In Figure (jl]) we show the temperature distribution in the equatorial plane for a simula- 
tion rotating with a period of 3— days. The sub-solar point on the day-side is characterized 
by an radially-extended nearly isothermal region with an effective day-side temperature 
T d ~ 1200A'. Upon adjusting to a hydrostatic equilibrium, a slightly negative temperature 
gradient is established so that the reprocessed stellar radiation can penetrate into the plane- 
tary envelope. Nevertheless the day-side photosphere is essentially isothermal with a density 
profile 



where is a planet's radius at the base of the isothermal region and g = GM v jr\ is the 
surface gravity of the planet. The imposed stellar heating falls off as a function of longitude, 
thus the radial extent of this isothermal region decreases with the inclination angle between 
the local zenith and the position of the host star overhead. 



Below the isothermal photosphere, a cooler region at lower depths is also evident on the 
day-side, associated in part with a cool return flow from the night-side. Figure ([5]) shows 
temperature distribution both throughout the entire photosphere, and focusing on structure 
on the night-side. Despite the added effect of rotation, a clear day-night delineation is 
still apparent, with the night-side characterized by effective temperatures T n from 310 to 
500K. The average night-side temperature is 380K, slightly smaller then the non-rotating 
simulation with the same opacities. This slight decrease in average temperature is due 
to increased cooling associated with rotationally modified flows discussed in the next sub- 
section. Slightly hotter regions near the terminators, associated with jets from the day-side, 
are apparent with temperatures reaching ~ 500K. 



3.2.1. The Day-Side Isothermal Surface 




(7) 



3.2.2. Azimuthal Effective- Temperature Distribution 
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Fig. 4. — The temperature distribution at the equator of a planet rotating with a period 
of 3days. The day-side is characterized by a large isothermal area near the top, the extent 
of which falls off with increasing longitude. A cool region at depth is also evident due to 
the combined effects of a negative temperature gradient from the internal heating, a positive 
temperature gradient from the irradiation, and a cooling return flow from the night-side. 
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Fig. 5. — The temperature at the photosphere of a planet rotating with a period of 3days. 
The upper panel shows the distribution over the entire planet, while the lower panel highlights 
the temperature structure on the night-side from = | to = A clear day-night 
delineation persists, despite complicated dynamical structure, due to substantial radiation 
near the terminators. 
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3.2.3. Thermal Current and the Corilois Effect 



The upper panel of Figure ([6]) shows the velocity magnitude at the photosphere (\v\ 



v \ + v e)- As i n the non- rotating model, material is moving quite rapidly, reaching speeds 
of ~ 4km/s near the terminators. Eastward (prograde or > 0, i.e. in the same direction 
as the unperturbed spin) moving material appears to be funneled from the sub-tropical 
latitudes (16*1 > 0) into an equatorial jet near = | , while westward (retrograde or < 0) 
moving material is pushed from the subtropical zone toward the poles near = It is 
these flow structures, rapidly advecting energy from the day-side, that account for the hotter 
regions seen on the night-side in Figure (jl]). To understand this structure, we must evaluate 
the 9 component of the Coriolis force, given by ^Qv^sin (8), which is shown in the lower 
panel of Figure Although material near the equator doesn't feel any Coriolis force in 
this direction, it is clear that material at higher and lower (sub-tropical) latitudes does. The 
asymmetry imposed by the rotation (i.e. the azimuthal velocities) causes the fluid moving in 
eastward and westward directions to behave significantly different then in the non-rotating 
case. 

The approximate magnitude of velocity can be estimated from equation ([2]), considering 
only the pressure gradient term. Assuming an approximately constant acceleration around 
to the night-side of the planet given by a ~ — ^VP, the velocity at the terminator should 
be given by 

"2 7 R G 



Vt 



1/2 

(T d -T n ' 



Given a day-side temperature of T d ~ 1200 AT and an average night-side temperature of 
T n ~ 350A', flows should achieve Mach numbers of ~ 2 near the terminator. The sound- 
speed ^RcT/yu) 1 / 2 ~ 1.7 km/ s at the terminator, yielding a local Mach number, as predicted, 
of ~ 2. 

Also evident in the plot of velocity magnitude in Figure is the marked decrease in 
velocity where the eastward and westward flows converge. Neither flow is able to instigate 
circumplanetary flow at the surface. Figure ([7]) shows an equatorial slice of the azimuthal 
velocity at the equator (upper panel) and at higher latitudes (lower panel). It is evident 
from this plot that the eastward moving flow does continue around the planet at depth near 
the equator, while the westward moving fluid continues around the planet at higher (and 
lower) latitudes. Because of the effects of rotation shown in Figure (Q, the convergence 
point is near 4> = for the equatorial flow, and near 4> = ^f- for flows at higher and lower 
latitude. This flow pattern implies that, upon converging, one of the two flows has undergone 
substantially more cooling then its counterpart. Thus, at the equator, when the eastward 
flow encounters the westward flow near = ^ (past "mid night"), the former is cooler and 
sinks below. The opposite is true at higher latitudes, with the westward flow experiencing 
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more cooling and sinking below the eastward flow. This cooling trend can be seen in the 
lower panel of Figure Lastly, very little motion is apparent deeper in the planet, as 
the flows are confined to a relatively small region near the top of the planet, supporting our 
assumption of a spherically symmetric inner boundary condition. 

3.2.4- Sub-surface Thermal Stratification 

Another method of visualizing the structure of the planet is through pressure-temperature 
profiles. Of crucial importance when calculating the emergent spectra, the pressure-temperatures 
at four different longitudes are shown in Figure (jHJ) for both the equator and 9 = 35°. All 
profiles agree well at pressures above O.lbar, again supporting the assumption of spherical 
symmetry at depth. However, below this pressure (or above this height), their behaviors are 
quite different. 

The day-side profile undergoes a significant temperature inversion (i.e. temperature 
begins to increase with radius and decrease with pressure). This temperature gradient allows 
radiative diffusion from the photosphere to the planetary interior. This excess flux is advected 
to the night-side deep down in the planetary envelope. Comparing to Figure it is evident 
that the lowest temperature region, ~ 650K, is associated in part with the cool, return 
flow. This transitional region is analogous to the thermocline in the terrestrial ocean which 
separates the surface and deep water layers. 

Near the upper atmosphere (where the pressure is low), the temperature distribution 
contains a large, approximately isothermal region on the day-side, clearly dominated by the 
stellar irradiation. In contrast, the temperature of the night-side (0 = ir) monotonically 
decreases throughout the entire atmosphere, reaching a temperature of ~ 300K at the pho- 
tosphere. While the day-side is fully radiative, there is a region near the photosphere on the 
night-side that is convectivly unstable in the radial direction. As a consequence of efficient 
convective transport, the P — T distribution on the night side is approximately adiabatic. 

The profiles near the terminators (at both = f and = ^) exhibit more complex 
structures due to considerable differences in the advective transport of heat. They also show 
differences between the equatorial values (left-hand panel) and higher latitudes (right-hand 
panel). At the equator, the temperature at = | exhibits an isothermal region ranging 
from 0.1 to 10~ 3 bars. This terminator is associated with the prograde flow. By the time 
this flow reaches = ^ it has cooled substantially. A slightly hotter region is evident from 
the westward flow at lower pressures (larger radius) near the photosphere. For the profiles 
from 9 = 35°, it is at = ^ where an approximately isothermal region exists at depth. 
The = f profile at high latitudes decreases monotonically, as very little heat is advected 
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Fig. 6. — The upper panel shows the magnitude of the velocity at the photosphere, given 
by \v\ = yjvl + vj. The lower panel shows the latitudinal component of the Coriolis effect 
(— 20v^,sin(^)) at the photosphere of a planet spinning at 3days. The direction of the 
latitudinal Coriolis force is different for eastward and westward moving material, causing 
eastward flows to be focused toward the equator, while westward flow is funneled toward the 
poles. 
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Fig. 7. — The velocity distribution at the equator (upper panel) and a latitude of 6 = 35° 
(lower panel) of a planet rotating with a period of 3days. At the equator, eastward flow 
(red) is able to circumnavigate the planet at depth, while at higher and lower latitudes, it 
is westward flow (blue) that is able to traverse all the way around the planet. The sinking 
of one flow under the other is due to different degrees of cooling due to rotationally altered 
dynamics. 
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eastward. 



3.3. Opacity Effects 

The effect of opacity can not be understated. Opacities regulate the efficiency of both 
the absorption of the incident stellar irradiation on the day-side and the re-radiation from the 
night-side. As noted by previous authors, major uncertainties in composition, metallicity, 
and chemistry all cause significant changes to the opacity. Furthermore, the inclusion of 
clouds, characterized by models of particle size distributions and vertical extent, tends to 
smooth out wavelength dependent opacities, resulting in spectral energy distributions that 
more clos ely approxima t e bla ckbodies. In our current models, opacit ies are found using the 



tables of iPollack et al.l (119851 ) for lower temperatures coupled with lAlexander fc Ferguson 



(119941 ) for higher temperatures. These are Rosseland mean opacities and include the effects 



of atomic, molecular, and solid particulate absorbers and scatters. 

In light of the uncertainties associated with opacity, we have studied the effect of uni- 
formly changing the opacity by some multiplicative factor. More detailed studies of specific, 
temperature and density dependent augmentations to the opacity will be presented elsewhere. 
Figure shows the temperature distribution across the photosphere of a planet with opac- 
ities reduced by both a factor of 100 (upper panel) and a factor of 1000 (lower panel). In 
comparison to Figure ([5]), night-side temperatures are significantly higher for both models 
with lower opacities, and distributions are smoother. In addition, due to changes in the flow 
detailed below, the hottest spot is displaced slightly from the sub-solar point. At the equa- 
tor, the displacement is ~ 10° and ~ 20° eastward (in the direction of rotation) for opacity 
reductions of 100 and 1000 respectivly. This displacement is largest at the equator, with 
maximum temperatures at higher/lower latitudes occurring closer to the sub-solar longitude. 



In a previous analysis (IBurkert et al.l (120051 )). we derived a formula for the night-side 



temperature by equating the radiative timescale with the crossing timescale. This two-zone 
(day-night) model assumes that: the advective heat flux is larger then the heat flux from 
the interior, the heat carried by a day-night thermal current is determined by the amount 
of radiative diffusion during the hemispheric circulation, and the night-side radiates all the 
heat advected to its proximity as a black body. With these assumptions, the night-side 
temperature can be estimated by 



\3irK d aR p J 



1/4 

(9) 



which is an decreasing function of opacity. In this formula, v is the average advective speed 
which is on the order of the day-side sound-speed c<j, and n d is the day-side opacity. The 



18 



T(K) 



TOO 



Fig. 8. — The pressure temperature profiles within a planet spinning with a 3day period. 
The left-hand panel shows profiles at the equator (6 = 0), while the right-hand panel shows 
profiles from a latitude of 9 = 35°. The individual lines represent = (solid line), <p = | 
(dotted), = Ti (dashed), and (f> = (dash-dot). The dot on each profile denotes the 
location of the photosphere. 
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Fig. 9. — The temperature at the photosphere of a model with the opacities reduced by a 
factor of 100 (upper panel) and 1000 (lower panel). Lower opacity fluid absorbs the incident 
stellar irradiation deeper in the atmospheres. This higher density material is able to advect 
the energy to the night-side more efficiently, leading to larger night-side temperatures. 
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increase of night-side temperature with decreasing opacity reflects the depth that incident 
stellar irradiation is deposited on the day-side. If the atmosphere contains grains with an 
abundance and size distribution comparable to that of the interstellar medium, only shallow 
heating occurs on the day-side, and the circulation does not effectively transmit heat to 
the night-side, which then cools well below the day-side. However, as the abundance of 
grains in the atmosphere is reduced, the stellar radiative flux penetrates more deeply into 
the atmosphere on the day-side, and the higher density atmospheric circulation carries a 
larger flux of heat over to the night-side. 

For the parameters used here, equation ([9]) predicts T n ~ 250K for the interstellar 
opacity simulation, while for the lower opacity simulations, the formula predicts ~ 750K 
and ~ 950K for reductions of 100 and 1000 respectively. Inspection of Figures (J5J) and 
Figure (jUJ) show average night-side temperatures of ~ 380K, ~ 700K, and ~ 890K for the 
same three cases. These results clearly indicate that the night-side temperature decreases 
with the magnitude of the opacity. The predicted values from equation (Q generally agree 
with those from the simulation and it provides a framework for understanding the global 
heat flow. The differences in these results can be attributed to several factors; the surface 
heat flux carried by the thermal current from the day to night-side is not entirely radiated on 
the night-side, but rather cools as it travels and advection at depth plays an important role 
in transporting heat. These effects can be incorporated into a more comprehensive four-zone 
(day-night and interior-photosphere) model where we examine the energy transfer within the 
optically thick regions below the planetary photosphere. 

Before the presentation of the four-zone model, it is useful to analyze the results of the 
numerical simulation. These calculations indicate that changing the opacity of the atmo- 
sphere not only alters the night-side temperature, but it also modifies both the flow dynamics 
and interior structure. For large opacity (our standard case), the isothermal region on the 
day-side is relatively shallow, and thus the increased cooling leads to lower night-side temper- 
atures. The large temperature differential promotes a fast flow velocities around the planet. 
As the opacities are decreased, and the night-side temperature increases and the velocities 
decrease. In the lowest opacity simulation, the velocity remains subsonic throughout the 
entire simulation. A more uniform temperature across the planet surface also allows for 
circumplanetary flow near the equator by reducing the pressure gradient. Figure (fTUl) shows 
the velocity at the photosphere of the two low-opacity simulations. These results should be 
compared to the standard- value results shown in Figure As in the standard case, flow at 
higher and lower latitudes still travels westward at the ^ terminator and circumplanetary 
flow at the surface is suppressed due to increased cooling times. 

In Figure (jTTI) we show the pressure-temperature profiles at the equator for the two 
reduced opacity simulations. In contrast to Figure (jHJ) where the night-side was fully con- 
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vective, both lower opacity simulations exhibit isothermal regions in the upper atmospheres 
around the entire planet. Below the photosphere, there is a slightly negative temperature 
gradient so that the day-night advective heat flux deep beneath the photosphere can radia- 
tively diffuse to the planet's photosphere. As it is to be expected, the radial extent of the 
nearly isothermal atmosphere is largest with the lowest opacity. Both simulations retain a 
convective regions below the isothermal regions on the night-side. 

3.3.1. A Four-Zone Model for the Sub-Surface Thermal Structure of Planetary 

Atmospheres 

Discussions in the previous section clearly indicate that thermal currents not only trans- 
port heat from the day to night, but also from the night to day-side at other radii. In addition, 
the day-side atmosphere is radiative and isothermal, while the night-side is convective and 
adiabatic. Finally, quasi hydrostatic equilibria is maintained in most regions except the 
subduction zones which separate the thermal currents. We now take these effects into con- 
sideration in the determination of the night-side temperature and the sub-surface thermal 
structure. 

The two-zone approximation in equation Q is based on the assumption that the radia- 
tive and crossing time-scales in the planetary atmosphere are nearly equal. We show in the 
next section (Figure ffl2l) ) that this assumption is only marginally satisfied on the night-side 
near the photosphere. In the standard opacity model, the radiative timescale around the 
terminators is several orders of magnitude shorter than the crossing timescale. This allows 
for significant cooling, yielding much lower night-side temperatures then would be expected 
from equation (Q. Similar behavior is seen in the lower opacity models. 

In order to analytically account for the transport of heat well below the photosphere on 
either side of the planet, we divide it into four zones, representing the day- night and interior- 
photosphere regions. There are then 4 sets of thermodynamic variables including: Pds, Pdi, 
Pns, Pni, Tds, Tdi, Tns, and T^i where P and T are pressure and temperature, and the 
subscripts D, N, S, and I represent day, night, photospheric surface and interior respectively. 
There are three additional quantities which connect these regions: the velocity between the 
day and night interiors, v ac [ v , the radial distance between the photospheric surface and the 
planetary interior, Id, and the density of the gas at the planetary interior on the day-side 
Pdi- These quantities can be solved simultaneously with the following 11 equations. 

Although the total stellar incident flux F = L^/A-kD 2 (where is the stellar luminosity 
and D is the distance between the star and the planet), irradiates only on the day-side, the 
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Fig. 10. — Velocity at the photosphere of models with opacities reduced by a factor of 
100 (upper panel) and 1000 (lower panel). Velocities, determined by both the day-night 
temperature differential and the cooling efficiency, decrease with deceasing opacity. 
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Fig. 11. — The pressure temperature profiles at the equator of a planet spinning with a 3day 
period. The left-hand panel shows profiles for the run with opacities reduced by a factor 
of 100, while the right-hand panel shows profiles for a simulation with opacities reduced 
by a factor of 1000. The individual lines represent = (solid line), = | (dotted), 
<p = 7r (dashed), and = ^ (dash-dot). The dot on each profile denotes the location of the 
photosphere, considerably deeper then in the standard opacity simulation. 
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condition for thermal equilibriums implies that 



2a(T* s + T* s ) = F. (10) 

Hydrostatic equilibrium at the photosphere of both the day and night sides gives: 

2 9 n 2 g 



DS 



NS 



3 1 - 1VJ Q 

Kds <j i^ns 

For the simplicity of analytic approximation, we represent opacity as 

K = K T P . 



11) 



(12) 



In the temperature and density ranges which are relevant to the present investigation, the 
opacity is roughly independent of density, thus we adopt the approximation kq = 0.0391 and 
P = 0.641. 

Our numerical results indicate that the pressure may change by more then two orders of 
magnitude between the planetary surface and the temperature inversion layer. In contrast, 
the temperature changes by less than a factor of two. In the spirit of analytic simplicity, 
we adopt an isothermal approximation for the hydrostatic envelope when we determine the 
pressure at the planetary interior on the day-side, such that 



Pdi = Pds^xp 



kg 



(13) 



where c Sj d = (i^-gTds / f^) 1 ^ 2 is the sound speed. However, to calculate the temperature 
variation, we adopt a radiative diffusion approximation in the radial direction. Including the 
azimuthal advective transport, we have 



T, 



Dl 



DS 



4-/3 



(4-/3) 



PDlV a dvC sd 

F 



exp 



gk 



rpA 

1 NS 
1 DS 



where 



(14) 



(15) 



Pdi — P-Pdi/RgTdi- 

We choose Id to be the depth where advection carries half of the incident flux, i.e. where 

PDlVadvCld ~ F/2. (16) 

The advective velocity at depth is driven by the pressure differential, such that 
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(17) 
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We assume the consequence of the advective transport is to thoroughly mix the gas so that 
planetary interior becomes isothermal with 



T m = T DI (18) 

For the standard opacity model, the night-side remains fully convective. Assuming 
convection is efficient, the envelope of the night-side is expected to be adiabatic, such that 

7-1 

Tm = f Pm \ - (ig) 
Tns \ Pns J 

The energy equation for the night-side can be written as 

1 PDl\ _ Vconvg^Rp ^ 

Pni J RgTni 

where the convective velocity v conv can be obtained from the mixing length theory. For the 
low opacity models, the upper regions of the night-side become stabilized against convection. 
We can then replace the energy equation with the diffusion approximation and recalculate 
the thermal structure accordingly. 

With these equations, we can provide a set of algebraic equations which essentially 
reproduce the behavior of our numerical simulation. These equations are also more compre- 
hensive than that in equation (jUj). For opacity similar to that of the interstellar medium (i.e. 
no reduction in re ), the night-side remains convective so that equations (|T0|) -( l20|) are valid. 
In this limit, T^s is sufficiently larger than Tns such that equations (flOl) and (TTTT) lead to 
complete information on the surface layer of the day-side, 

T* s = F/2a, P DS = 2g/3K Tg s . (21) 

In comparison with the large variations in equation ffl~3l) . the day-side is nearly isothermal 
with Tdi ~ Tds- At the depth where half of the incident stellar flux is advected from the 
day to night-side, equation fTl6l) yields 

Pmv adv * f , l d c ^ln ( ^M&A . (22) 

In the fully adiabatic night-side, we find from equations ffTTj) . (j!3j) . (fT7|) . and fTl9l) that 

Pns — Pds{Tni/Tns) > (23) 

and 
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From equations (jTTl) . (fT8|) . and (|20l) . we find 



exp 



2< d 



^adv 



In the limit that t> a( fo < c Sj d, equation (|25|) reduces to 

('OH I' 



vl/3 



(25) 



(26) 



The convective speed can be estimated from the mixing length theory which is mostly de- 
termined by T NI ~ T DS and These quantities have very little or no dependence on k 
so that v a d v does not change significantly as a function of k,q. In this limit, equation ([2^ 
becomes 
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(27) 



Thus, the four-zone model generates a similar T/vs dependence on k as equation ([9]). In 
comparison with numerical results, equation fl27|) does reasonably well for high opacity con- 
vective simulations, but expression would be improved by considering the possibility that 
the night-side may also become radiative in the limit of very low n Q . 



4. Analysis and Model Comparisons 



Burkert et al 



A number of groups, inc l uding Showman fc Guillot ( 20021). Cho et al. ( 2003b . 

f l2005h . lCooper fc Showman! J2OO5UCI10 et al.l J2006h . and lLangton fc LaughlirJ J2007h have 
carried out non-linear numerical simulations studying the dynamics of hot- Jupiter atmo- 
spheres. Both the methods and results vary considerably. In this section, we attempt to 
compare their methodology to that presented here, concentrating on the assumptions that 
lead to differing results. 



Both IShowman &: Guillotl (120021 ) and ICooper &: Showman! (120051) solve the primitive 
equations, the former using the EPIC code developed by lDowling et al. (|l998l) and the late r 
using the ARIES/GEOS Dynamical Core model initially developed by ISuarez et al.l (119951 ). 
The primitive equations, cast in isobaric coordinates, are widely used in meteorology. When 
deriving them from the full Navier-Stokes equation, a key assumption of hydrostatic equi- 
librium is invoked. Assumption of hydrostatic equilibrium explicitly links the thickness of a 
layer to its local temperature. Radial motion can only occur in conjunction with divergence 
along isobars, the magnitude determined by that which maintains the condition of hydro- 
static equilibrium. With this assumption, the continuity equation and hydrostatic condition 
take the place of the radial momentum equation, and can be written as 

8U 0, (28) 



V p -v 



dp 
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where v is the horizontal velocity and u is the vertical velocity with respect to pressure 
coordinates. The resulting radial velocity is slow in comparison to the horizontal (or isobaric) 
motion. Although an excellent assumption for thin terrestrial atmospheres, it neglects critical 
flows present in thicker atmospheres such as those of hot- Jupiter's. Although deviation from 
hydrostatic equilibrium provides the dominate radial acceleration, high velocity azimuthal 
flows on the night-side give rise to non-negligible radial forces neglected in the primitive 
equations that contribute to the ability of the converging flows to pass under one another as 
seen in Figure©. 

Our numerical results and analytic approximation clearly indicate the dependence of 
the night-side temperature on the efficiency of radiation transfer. The heat flux carried by 
the thermal current is determined by the penetration depth of the incident stellar radiation. 
Radiative diffusion and convection are essential processes which regulate the heat diffusion 
well below the photo s phere of th e planet. In conjunct i on wit h the primitive equations, both 
Showman Sz Guillotl (120021 ) and ICooper fc Showman! (120051 ) utilize a Newtonian radiative 
scheme to approximate stellar irradiation. In principle, this scheme is only appropriate for 
energy deposition and emission in optically thick regions. The Newtonian radiative scheme 
provides a heating/ cooling term to the energy equation that relaxes the temperature toward 
some pre-defi ned equilibrium distribution on some radiative timescale. The forcing term can 
be expressed ( ICooper fc Showmanll2005l ) as, 



T(9,^p,t)-T eq (0,0, p) 

Teq (P) 



(29) 



For the radiative times cale. r fff . IShowman fc Guillotl (120021 ) assume a cons tant valu e given 



by that at 5bars, while ICooper fc Showman! (120051 ) use the calculations of llro et al.l (120051 ) 
to set a radiative relaxation timescale that is dependent on the local pressure. As noted by 
the authors, this approximation is crude, but allows rapid computation of a large number 
of models. In order for the Newtonian approximation to be viable, the radiative timescale 
{Trad ~ ^p-) must be much longer than the crossing timescale {r x ~ f^)- m this limit, 
the temperature distribution will be determined by the dynamics, rather than the assumed 
equilibrium distribution. This assumption becomes problematic in the upper atmosphere. 
In Figure (JT21 we plot the ratio of timescales Ir f A at the photosphere of the planet. 

It is clear from Figure (TT21) that T ra d < r x throughout the photosphere. The most 
significant deviation occurs on the day-side, where there is very little motion, and T raa - exceeds 
r x at the planetary photosphere by many orders of magnitude. Near the terminator, where 
the stellar irradiation falls drastically, the surface radiative timescale is also much shorter 
then the dynamical timescale; it is here that the deviations from a Newtonian radiative 
scheme are most important. The flow radiates a significant portion of its energy before it 
flows to the night-side, allowing for the sharp edges in the temperature distribution seen 
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Fig. 12. — The ratio of radiative timescale (r ra d) to the crossing timescale (r x ) at the photo- 
sphere of the simulation with interstellar opacities. Near the terminator, the primary region 
of interest, the radiative timescale is ~ 5 orders of magnitude shorter then the crossing 
timescale. This rapid cooling accounts for the sharp temperature gradient near the termi- 
nator seen in Figure (j5J). 
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in Figure ([5]), despite high velocity flows. In contrast, the Newtonian assumption that 
Trad >> t x , implies that the fluid carries a significant quantity of heat, leading to an overall 
distribution that will be more uniform with longitude. Finally, despite substantial radiation 
near the terminators, the near equality of T raa < and r x at <fi = it demonstrates that the 
temperature at the back-side is determined purely by the advection of energy by the flow. 

To again highlight the effect of opacity, Figure (TIB"]) shows the ratio of radiative and 
crossing timescales, for the runs with varying opacity. The upper and lower panels show 
the simulations with opacities reduced by a factors of 100 and 1000 respectively. Again, the 
most important area in determining the behavior is near the terminators. In comparison 
to Figure (TT2"j) . it is obvious that lower opacity fluid is able to advect energy to the night- 
side more efficiently due to increased cooling times. As mentioned above, lower opacity 
on the day-side, allows heat to penetrate further into the planet, and re-radiation near the 
terminator is somewhat suppressed. Although Figure ([TBI shows a marked increase in ^ L 
from the interstellar opacity simulation, the Newtonian approximation {r ra( i » r x ) is never 
satisfied. The largest values, ~ 3, only occur on the night-side in the high velocity 
circumplanetry jet. 



The studies of lCho et all (120031 ) , Menou et al.l (120031 ) , ICho et all (120061 ) , and lLangton &: Laughlin 

( 120071 ) take a d ifferent approach, concentrating on the effect of eddies and waves on the over- 
all dynamics. ICho et al.l (120031 ) iMenou et al. (120031) and lLangton &: Laughlin! (120071 ) use 



the shallow- water equations, while ICho et al.l (120061 ) solves the more generalized equivalent- 
barotropic equations. The equivalent-barotropic equations can be derived by vertically inte- 
grating the primitive equations described above. The resulting equations describe fluid flow 
in a single, isentropic layer whose scale height can vary. Stellar heating is prescribed by a 
deflection of the scale-height as a function of position. Concentrating on a single layer allows 
for the high resolution necessary to study the effect of turbulent eddies and waves. However, 
allowing for only a single layer implies that the atmosphere is radial isothermal. Although 
not a bad assumption in the upper regions on the day-side, the rest of the atmosphere ex- 
hibits significant radial temperature structure. As noted above, cooling in the upper regions 
of the night-side allows material to sink radial and initiate a return flow. Disallowing this 
flow would significantly alter the resulting dynamics. 



The final dynamical study of hot- Jupiter's was done by iBurkert et al.l (120051 ). The 
numerical model they used is quite similar to the model presented in this paper. They solve 
the full hydrodynamical equations, given by equations ((H), (J2J) and (jHJ) together with flux- 
limited radiation diffusion. However, they restrict their attention to the r — plane from 
= — ► 7r, neglect the curvature terms in equation (f 2 ]), and ignore the effects of rotation. 



Comparing our non-rotating results to IBurkert et al. 



(120051 ) 'Case 1' with our 'standard' 



opacity, we find that the two models agree quite well. The only substantial differences 
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Fig. 13. — The ratio at the photosphere of models with opacities reduced by a factor 
of 100 (upper panel) and 1000 (lower panel). Although lowering the opacity decreases the 
ability of the fluid to radiate its thermal energy, the radiative timescale remains several 
orders of magnitude shorter then the crossing timescale for most of the photosphere. 
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are that our backside temperature is slightly lower (~ 415K compared to their 480K), and 
the convergence point is seen to oscillate in the simulations presented here. Our cooler 
temperatur es can easily be expla ined by considering the increased compressional heating that 



the fluid in iBurkert et al.l (120051 ) simulations experiences as it hits the boundary at = 7r, 
and our extension into three- dimensions. Flows that are able to spread in latitude will cool 
more then those confined to the equator. Given that our non-rotating simulation results in 
two symmetric convective cells suggests that by restricting attention to <p = — > 7r, they 
did not miss any fundamental physics for that scenario. However, the addition of rotation 
significantly changes flow patterns and allows for increased cooling. 

In summary, four fundamental differences are included in the models presented here, 
that are not included in some way in the previous methods: self-consistent radiative trans- 
fer, solution of the full radial momentum equation, 3-dimensions, and rotation. Noticeable 
changes occur due to the inclusion of each. Self-consistent radiative transfer allows significant 
cooling as the fluid travels around to the terminators and night-side of the planet. Coriolis 
forces alter the structure of the flows, alternately compressing and diverging streamlines 
moving around the planet in different directions, and the multi-dimensional aspect changes 
the resulting temperature distribution. Finally, the full treatment of the radial momentum 
equation yields significant, non-hydrostatic radial structure and motion. Disallowing effi- 
cient cooling, and limiting radial motion, while subjecting the planet to the continual stellar 
energy input, will lead to much more uniform temperature distributions across the entire 
planet then those that we observe here. 



5. Discussion 



In this paper we have considered the flow dynamics and heat redistribution in a tidally 
locked hot- Jupiter orbiting a solar-type star with an orbital period of 3days. We utilized 
a three-dimensional model that solves the full hydrodynamical equations and models radia- 
tive transfer the flux-limited radiation diffusion. We show that the temperature distribution 
across the planetary atmosphere is a sensitive function of its opacity. Our models exhibit 
strong day-night temperature contrasts, despite strong winds, the size of which increases with 
increasing atmospheric opacity. Large temperature differences are due, in a large part, to sig- 
nificant cooling of the flow near the terminators of the planet. Rotational effect significantly 
alter the flow pattern, allowing for increased cooling and suppressing surface circumplanetary 
flow in the higher opacity atmospheres. 



In our standard model, opacitie s are calculated using the tabl es of iPollack et al.l (119851 ) 
for lower temperatures coupled with lAlexander fc Ferguson! (119941 ) for higher temperatures. 
These are interstellar Rosseland mean opacities and include the effects of atomic, molecular, 
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and solid particulate absorbers and scatters. Although there is significant room for im- 
provement in our treatment of opacity, we explore the effect of uniformly red ucing them by 



factor s of 100 and 1000. In agreement with the two-dimensional simulations of lBurkert et al. 



(120051 ). we find that lower opacities yield higher night-side temperatures. Lower opacities 
allow incident stellar irradiation to be deposited at larger depths on the day-side, decreasing 
the cooling rate as this energy is advected to the night-side, and increasing night-side tem- 
peratures. The relation between the night-side effective temperature and the opacity in the 
planetary atmosphere is verified by a more comprehensive four-zone analytic model which 
takes detailed account of the effect of radiation transfer on both the day and night-side of the 
planet. This model highlights the importance of appropriately treating radiation transfer in 
the simulation of atmospheric dynamics on close-in gas giant planets. 

The actual value of opacity in the atmospheres of these planets is highly uncertain. 
However, the day-side temperature of short period gas giant planets around solar type stars 
is Tds ~ 1200K, which is below the grain destruction temperature. We note that the flow 
returning to the day-side is dredged up from a cooler interior {Tui < Tds) at the sub- 
solar point. Thus, throughout the thermal circulation, refractory magnesium-rich grains are 
preserved. Other species of less refractory silicate grains will sublimate near the planetary 
photosphere on the day-side, condense on night-side surface, and be carried along with the 
returning current beneath the surface. These thermal currents become turbulent as they 
generate specific vorticity and excite instabilities. Small grains are well coupled to the gas 
and turbulence induces them to collide and coagulate. In a static atmosphere, particles with 
sizes s p and density p p attain attain terminal velocities of 



v t ~ {gp P Sp/pg) 1/2 
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in a background gas with density p g . Since the distance over which these particle attain their 
terminal speed, L term ~ P P s p /p g , is smaller than the density scale height ~ c 2 s d /g, particles 
larger than ~ 1mm cannot be carried by the upwelling current on the day-side, and are 
left well beneath the planet's photosphere. This potential channel for grain-gas separation 
implies that it is very likely that the opacity in the planet's atmosphere is well below that 
of the interstellar medium. A detailed analysis of the grain evolution in the atmosphere of 
short-period gas giants will be presented elsewhere. 

We also presented a det a iled compa r ison o f our model to previous dynamical models 



of Showman fc Gui 



lot 



l2002iuCho et al.1 (I2003h. Burkert et all ( I2005f ). ICooper fc Showman 



( 120051 ). ICho et al.l (120061 ). and lLangton fc Laughlinl (120071 ). The approaches among these 



groups varies greatly, and we have attempted to highlight several of the crucial differences. 
The most significant differences that we include are solving the full fluid equations in all 
three-dimensions, rotation, and our treatment of radiation via the flux-limited diffusion 
method. Flux-limited radiation transport allowed for a self-consistent treatment of the flow 
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of radiation throughout the planet. This treatment produced significantly shorter radiative 
timescales then previously calculated, yielding much lower night-side temperatures. 



The dynamics and heat distribution within the atmospheres giant planets allow us 
to probe fundamental questions surrounding planet formation. The diversity of planetary 
sizes must ultimately result from interior structure variations arising from formation and 
evolutionary processes. These atmospheres not only serve as a valuable observable links to 
th e interior, but may al so pl ay a role in regu l ating p lanetary contraction rates. As suggested 
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Burrows et al.l (120001 ) and iBurrows et al.l (]2006al ) an isothermal atmosphere will result in 



decreased heat transfer and increased planetary contraction timescales. However, a strong 
day-night temperature difference may allow an avenue for cooling via the night-side. Based on 
the simulations presented here, it is clear that the efficiency of redistribution decreases with 
increasing opacity. High opacity atmospheres, similar to interstellar values, will have cool 
night-side temperatures, providing an avenue for internal heat loss and radial contraction. If, 
as is expected, there is a significant amount of grain growth and sedimentation, thus reducing 
the atmospheric opacity, the planets night-side will also contain a large isothermal region, 
possibly allowing for the retention of internal energy. Diversity in atmospheric opacity may 
lead to diversity in planetary radii. With the promise of improved observational techniques, 
including transit spectroscopy and full phase light-curves, the relation between opacity and 
temperature differential should be testable in the near future. 
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